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Abstract 



We re-examine the systematic errors in determination of the primordial 
helium abundance, Yp. We find that the systematics are significantly larger 
than the statistical errors. The uncertainty in (the determination of) Yp, is 
thus, larger than is currently claimed. Furthermore, most of the systemat- 
ics lead to underestimate of Yp . The new upper bound allows cosmological 
models with no non-baryonic dark matter in which ^baryons = ^BBN = ^dyn- 

I. INTRODUCTION 

The current limit on ftsBN^ the baryonic density parameter that matches big bang 
nucleosynthesis, is flsBN < 0.06 h'^Q [1] (Where h^o is the Hubble constant in units of 50 
km s~^ Mpc~^.) This limit is in conflict with other observations of the density parameter 
0; From inferred motions of virialized systems on small scales we get a dynamical estimate 
to which is ftdyn = 0.1 — 0.2 [2]. In order to satisfy the above BBN limit one has to 
invoke unknown non-baryonic dark matter. Theoretical prejudice, based on the inflationary 
scenario, favor a flat universe with ftfotai = 1- However, the observed ratio of matter+gas 
to total mass in clusters [3] ^baryons I ^totai ^ 0.01 + O.ld/i^g' is in conflict with the above 
constraint on ftsBN for ftfotai = 1- The primeval baryon isocurvature model (PBI) [4], with 
^totai being only baryonic between 0.1 — 0.2, i.e. ftfotai = ^bbn^ seems to fit observations 
better than any other cosmological models [5], except for one - the observation limit on 
^BBN- This uncomfortable low limit has provoked many attempts [6] to increase the bound 
on ftsBN- Most of these trials have been unsuccessful [7], while, the few successful scenarios 
[8] [9] are based on drastic modifications of the standard picture. 
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In this paper we reexamine the observational basis for the constraints imposed on ftsBN- 
Of all the BBN isotopes, ^Li is the easiest abundance to determine, however, it is very fragile 
which makes it hard to relate it to its primordial content. Just recently the ^Li abundance 
was increased by a factor of 7 [10] as a result of studies of both young and old population 
stars which require the inclusion of diffusion as well as rotational mixing in the stellar models 
[11]. This increase can be constrained by further detections of ^Li in old stars, but progress 
there is still slow. Independently, improvements in reaction rates and the neutron half-life 
[12] [1], in addition to recent correction in the Teff-scale of halo stars [13] brought the ^Li 
abundance even further up. The new values for ^Li require higher values for the current "^He 
abundance for consistency. 

Deuterium, can only undergo net destruction following the epoch of BBN. Because of its 
steep dependence on the ratio of baryons to photons, ?y, the D abundance could play a major 
role in studying BBN, if constrained more strongly. Unfortunately, current values of Dp from 
interstellar medium determinations have a large scatter around D/H sa 1.5x10"^ [14], while 
the Solar System measurements are limited by the need to know the solar "^He abundance 
[15] which has not been derived accurately yet. When the uncertainties are accounted for, 
the current lower limits on D and D-|-"^IIe allow our upper bound on the helium primordial 
abundance of yp<0.255; however, a Ip higher than that runs into conflict with the presently 
estimated values of Dp [16]. Quite recently two groups reported the possible detection of 
extragalactic D in the spectrum of a high-redshift QSO [17] [18]. They derive a deuterium 
abundance of D/H sa 25x10"^ corresponding to Yp near 0.235, but caution that the observed 
feature has a high probability of being due to hydrogen absorption along the same line of 
sight. However, if the absorption is due to deuterium the derived abundance is so high, that 
it is in conflict with the solar system observations of D and "^He [19]. 

As a result of the recent uncertainties in the determination of ^Li and D, the abundance 
of "^He continues to play a key role in defining the observational constraints on ftsBN- The 
uncertainties in the BBN observational limits are dominated by systematic errors which 
generally escape the attention of those who use Yp limits for theoretical constraints. During 
the last few years there has been substantial improvement in the atomic data for helium, as 
well as of the methods for solving detailed radiative transfer. In this paper we take advantage 
of both in order to evaluate the systematic errors in the derivation of the primordial helium 
abundances. Our goal here is to obtain a quantitative estimate of the systematic errors and 
new upper bounds on Yp and on ftsBN^ respectively. We expand on the implication of this 
new bound on various cosmological scenarios. 



II. THE ^HE LIMIT 

The observational determination of Yp relies exclusively on the use of the optical recom- 
bination spectra of metal-poor, extragalactic ionized nebula - H II regions and H II galaxies 
[20] [21] [22]. The major advantage of using such H II regions is the theoretical simplicity of 
recombination radiation. In the absence of complicating effects, the He/H abundance ratio 
could be deduced from relative intensities of hydrogen and helium lines using theoretical 
effective recombination coefficients. The procedure of deriving the abundances from the ob- 
servations consists of adopting a simple Case B radiative recombination theory [23] followed 



by a set of "corrections" which are apphed to counter its inadequacy [21]. In Case B one 
assumes that the hne photons (mostly resonance lines of abundant ions) are scattered so 
many times that their downward radiative transitions can effectively be omitted from con- 
sideration. The other extreme assupmtion, Case A, which is valid for most of the emission 
lines observed in nebulae, is that all emitted line photons escape and there is no radiative 
transfer problem. The real situation, however, is intermediate (Osterbrock [24], Ch.4), non- 
linear and dependents on the ions involved. The important sources of errors inherent to 
the H II galaxy method have been already identified and discussed (see e.g., references [25] 
[26] and [21]). Most of these are of systematic nature and fall generally into the category of 
radiative transfer problems. Below is a list of some of these radiative transfer problems: 

(1) The invalidity of Case B for H . 

(2) The invalidity of Case B for He I. 

(3) The amount of neutral helium. 

All the listed effects, plus some we do not list (e.g., consideration of underlying stellar 
absorption), always lead to a systematic underestimate of Ip (see Skillman & Kennicutt 
[25]). 

We examine the above systematic effects with the help of detailed radiative transfer 
calculations for H, He I & II. We solve simultaneously the coupled statistical equilibrium 
and transfer equations in plane-parallel illuminated slab geometry using the approximate 
A-operator method of Scharmer [27]. This is a full multi-level non-LTE calculation for the 
atomic models of H, He I, and He II. We take special care of the helium resonance lines 
transfer, e.g.. He II La 304 A, which can be an important ionization source. In constructing 
the slabs we used photoionization equilibrium calculations with a mean escape probability 
formulation after Proga et al. [28] and similar to model ^8 of Clegg & Harrington [29]. 
The helium atomic model is similar to that of Almog & Netzer [30], but contains some 
new atomic data from the Opacity Project database, the new term values from Drake [31], 
and the coUisional rates evaluated by P.J. Storey from the cross sections of Berrington & 
Kingston [32] (see Sasselov & Lester [33] for some details on the atomic data). We compare 
the results of our computation to the calculations using Case B analysis with the standard 
"corrections" for the same physical problem. The comparison gives us a realistic estimate 
of the systematics in the determination of Ip from the observed data. 

We apply the results of our calculations of systematic effects to four recent values of 
Y and Ip as derived by different observers: Ip by Pagel et al. [21] derived from He vs. 
and He vs. N; Ip by Olive & Steigman [34] from a new set of low-metalicity H II regions 
observed and analyzed by Skillman et al. [22] and combined with the set of Pagel et al.(for 
both and N); and the abundance determination in the most low-metalicity system - I Zw 
18 by Skillman & Kennicutt [25]. 

We proceed by describing the specific problems and their overall contribution to the 
individual systematic errors for each case. 

[(1)J The invalidity of Case B for H : Several processes complicate the radiative transfer 
of hydrogen in H II regions [35], [36], [37]. First, the low-metalicity H II regions have high 
enough electron temperatures that some of the observed Ha fiux may be due to coUisional 
excitation. The rate depends on the ground-state population, N//o, which cannot be derived 



from recombination theory alone. Worse still, the amount of neutral hydrogen is unohservahle 
in H II regions, due to the lack of any spectral diagnostics. The ground-state population, 
N/jo, can be obtained from detailed radiative transfer calculations with geometry constrained 
by high-resolution observations. 

To demonstate the magnitude of this effect we consider NW H II region in I Zw 18 for 
which high resolution data is available. In this case we can calculate the equilibrium ratio 
of Nffo/Nff+ by considering the ionization balance from the known size, estimated Tg, and 
a Zanstra guess for the ionizing-photon supply rate [26]. The recent HST spectra at H Ly 
a of I Zw 18 [38], provides an additional constraint on the Lyman continuum flux of about 
2.9xI0^^s~^. We use our illuminated slabs results together with Tg and N^ from Skillman 
& Kennicutt [25] to get: Nffo/Nff+ sa 8 X I0~'^f2fi ■ Here /i and /2 are volume filling 
factors (from 0.5 to 0.0 1 - [39]), accounting for clumping and other non-uniformities in the 
gas density and optical thickness. We find that even such small amounts of i7° lead to a 1% 
to 2% effect on the derived helium abundance, via its effect on the emergent Ha fiux (see 
Table 3 in ref. [35]). 

Though small, the error is systematic; always resulting in an increased Y. The case of I 
Zw 18 is given above for illustration, but in general additional knowledge of the geometry of 
the H II region are necessary in order to build a realistic model [40], and such is impossible 
to obtain for distant extragalactic H II regions. 

A second process is incomplete absorption of hydrogen Lyman lines due to the finite 
optical thickness of the gas and dust, or photoionization of metals [36]. Balmer line emis- 
sivities are reduced relative to those for Case B. We calculate continuum and line opacities 
in order to use the results of Hummer & Storey [36] for Ne = 100 cm~"^. For I Zw 18 the 
effects due to dust are increasingly negligible due to its low dust-to-gas ratio [38]. Hence, 
we do not take them into account in our calculations of the systematics for this case. 

Another effect, which we do not include here, is due to dust emission. We point out, 
however, that for all the H lines of interest, it has a systematic deviation of the same sign as 
the effects already discussed. It is mostly significant for the radio recombination lines [36]. 

The overall systematic effect on the abundances based on Ha/(HeI A 6678) and Tg in 
the 16,000-20,000 K are in the range of I%-3%. 

[(2)]: The invalidity of Case B for He I : In general, recombination theory is rarely simply 
valid for neutral helium, because of the distinct singlet and triplet species of He I, with the 
strongly metastable 2"^S state of the triplets (Osterbrock [24], Ch.4.6). In the calculations of 
effective recombination coefficients the triplet series is often treated only as Case B, while 
the singlets are treated in both Case A and B (e.g. ref. [42] [43]). Therefore, whatever the 
detail and precision employed in such models, they do not account for important radiative 
transfer effects within He I, like the intersystem transitions coupling the singlet and triplet 
series (the latest models [43] include only those between n = 2 levels). Through this coupling 
the accurate calculation of the singlets of He I becomes dependent on the detailed treatment 
of the triplets, which are more sensitive to coUisional and self-absorption effects in low- 
density plasmas. The use of simple recombination theory excludes also any consideration of 
inter-species radiative transfer effects. One example is the coupling between the hydrogen 
Lyman continuum (A 912) and the He I (A 504) continuum, which can become optically 
thick only in accord with the optical depth of hydrogen, T912, because the dominant opacity 
at 504 A is due to the Lyman continuum. Once the Lyman continuum becomes optically 



thick, the ionization equihbrium of hydrogen shifts in favor of H I and affects Ne] hence — 
the total recombination rate of He I. 

In order to illustrate the problem and derive some estimates on the systematics involved, 
we have calculated detailed models, extending to lower densities the models of Almog & 
Netzer [30]. In Figure I our results are compared to Case B calculations, and to two 
observed low-metalicity H II galaxies. This line-ratio diagram juxtaposes the singlets (He I 
A 6678) to the triplets (He I A 7065). Through the fluorescent enhancement of A 7065, this 
diagram can also indicate large optical depth, T3889, in the He I A 3889 line. In the lowest 
density regime (I — 100 cm~"^), both the singlet A 6678 line and the He I A 5876 line are not 
sensitive to Te or A^e, and their ratio is practically constant. This is not true for the triplet 
A 7065 transition which is fluorescent (coupled to A 3889), its emissivity increases with T3889, 
as collisions and optical depth become more important. 

Figure I illustrates the degree to which existing models account for the observed fluxes 
of He I lines used in the abundance determinations. No model fits the observations well but 
the Case B models fail most conspicuously to come close to the generous observational error 
boxes. The singlet line He I A6678 is more reliable but still suffers from systematic effects. 
From Figure I it is also obvious that Case B models underestimate the lines sensitivity to 
Ne (for both He I systems), as is inherent to Case B. Notice how in the detailed models the 
density sensitivity increases with Te, as expected, while no such trend is discernible in the 
Case B models. 

Despite the availability of detailed H II region photoionization models [40] the current 
quoted values of Yp are based on assumption of Case B and the emissivities of Brocklehurst 
[44]. The latter neglect the metastability of the 2"^S state of He I, and are known to have 
other problems with the triplets fiuxes [42], thus affecting the results of analyses which used 
mostly triplet He I lines to determine Y. Unfortunately, although small, the discrepancy is 
systematic (underestimating Y). 

In summary, the systematic effect due to problem (2) - the invalidity of Case B for He I 
[44], would increase Y by 4%-6% (for He I A 6678) and 5%-8% (for He I A 4471) for T^ 
between 10,000 and 20,000 K. 

[(3)] Neutral helium: This problem arises from the lack of observational indication of 
the amounts of neutral helium (similarly to hydrogen) in an H II region. Yet any A^He" 
not accounted for underestimates the final abundance Y. Very detailed model calculations 
constrained by the best for any H II region observations, estimate A^He" at 3.2% in the 
Orion Nebula [45]. This nearby H II region is rich in metals and dust; dust extinction 
tends to keep ionization high and A^He" low. As we go to the extragalactic H II regions, the 
problem becomes extremely difficult to solve correctly and consistently, because of its strong 
dependence on geometry and composite H II regions [40]. The systematic error due to the 
presence of neutral helium can vary from zero to 3%. This follows from the range of values 
for the amount of neutral hydrogen (Problem I) we obtain from the available parameters of 
extragalactic H II regions. 

We sum up all the systematic effects, for each of the observations, taking into account the 
coUisional corrections that have been made already while deriving each of the abundances, 
i.e. we do not double count systematics. The helium abundance derivations in the last 
two studies [22] [25] [41] in Table I are based on (the more reliable) He I A 6678 line and 
Ha, hence we consider their systematics. For the study by Pagel et al., we have considered 



also the systematic effects on He I A 447f, A 5876, and H/3. The corrections for each Y 
value due to each radiative transfer problem are listed in Table f . The final upper bounds 
result from the systematic effects added, observers' coUisional corrections subtracted, and 
the 2(T observational statistical bound included. For each upper bound on Ip we derive a 
corresponding bound on the baryon-to-photon ratio, rj^edian in terms of rjio = ?yxfO^°, from 
the most recent theoretical relation of Kernan & Krauss [1] and their 95% confidence limits, 

V95%c.l.- 

Table 1: The Upper Bounds on Yp, rj in terms of rjio = ?yxlO^°, and ftsBN- 
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III. THE IMPLICATIONS FOR n 



From the above analysis we put most weight on two of the derived upper bounds: the 
one by Olive & Steigman, and the one of the NW H H region of I Zw 18 [25]. Their essence 
is different. While Olive & Steigman give a value of Yp based on careful extrapolation of the 
relation through the set of most recent H H region abundances, there is only a Y value for 
I Zw 18 which is by no means truly primordial Yp. However, K in I Zw 18 is very carefully 



derived from modern high-quality data by Skillman & Kennicutt [25], hence we are more 
confident about the size of the systematic errors (0.015) apphed to it. The gas in I Zw 18 
is not primordial, but is the closest to primordial we can get in H II regions [38], and it 
is good to keep its Y value in mind. As a result of all this, we suggest that current "^He 
abundance derivation methods do not exclude a Ip < 0.255, and perhaps even a value as 
high as Ip < 0.258, see Figure 2. We would like to point out that Ip < 0.255 is in agreement 
with the values and bounds for "^ Li and the current low limit on D and D -\-^ He (see section 

I)- 

We calculate an upper bound on ftsBN using the connection between rjio and ftbaryon- 
^baryonh"^ = 0.00365 TjiQ T2,^26 (here h is the Hubble constant in units of 100 and T2.726 the 
temperature of the CMBR in units of 2.726 K). From the above limits on Yp we get that 
flsBN can be as high as 0.22 h'^Q , and it might be even as large as 0.32 h'^Q . Our results 
using ?7io,95%c./. are summarized in Table 1. 

The new upper limits on ftsBN enable us to reconsider many cosmological models which 
were in confiict with the old limits. A universe with ftfotai = ^bbn = ^dyn is in agreement 
with the new bound and there is no need to invoke non-baryonic dark matter. One such 
a model is the PBI [4] model with sa 0.1. Moreover the model that was suggested by 
Gnedin and Ostriker in 1992 [8] in which ftsBN = ^dyn needed high values of Ip. Since 
this model was inconsistent with observations at that time, Gnedin and Ostriker had to 
turn to unconventional scenarios in order to justify the high Ip values. These Ip values are 
well inside our new limit and therefore this model can rely on conventional cosmology. And 
finally a value of Ip near 0.251 is all that is needed in order that the fraction of baryons 
in clusters [3] will agree with the value of ftsBN in a fiat universe dominated by cold dark 
matter. 



IV. CONCLUSIONS 

In this paper we do not determine a new value for the primordial helium abundance, Yp; 
we determine only the size of some important systematic errors. When accounted for, they 
provide an upper bound on Yp which is less confining than currently assumed. The new 
bound on Yp provides an upper limit to ftsBN which is much less restrictive for cosmological 
models. 

The new relatively high upper bound on ftsBN removes the need to turn to exotic models 
in order to reconcile between various observations of the density parameter. Moreover, the 
simplest and most esthetic model in which ftfotai = ^bbn = ^dyn ~ 0.1, a model that does 
not rely on unobserved entities (dark matter) or unconventional physics and is consistent 
with many observations ( [46] and references there), is also consistent with BBN limits. 

Because of the logarithmic dependence of Ip on ?y, a very accurate determination of Yp 
is needed in order to determine ftBBN precisely. Unfortunately, because of the observational 
limitations of the H II galaxy method, all we were able to calculate in this paper are upper 
bounds on Yp and, therefore, we are unable to constrain further the cosmological models. 
Two of the important sources of systematic errors in the H II galaxy method discussed 
here are related to the geometry and stellar content of the H II regions. Such constraints 
cannot be achieved with the resolving power of any current telescope. Fven with the best- 
observed nearest H II region (the Orion Nebula) the uncertainties in the helium abundance 
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are still difficult to reduce below 7% [45]. Therefore, modern sophisticated radiative transfer 
techniques cannot be applied with equal success to the extragalactic H II regions. As an 
alternative, extreme halo stars in our Galaxy may provide a way to derive Ip with precision 
interesting for cosmology. Much of the improvement will come from the fact that halo 
dwarfs are much closer to primordial abundances in their unprocessed atmospheres. This 
is illustrated in Figure 2, where some of the program halo stars are marked with respect to 
the H II regions of the lowest known metalicity. A project to derive helium abundances in 
extreme halo dwarfs and obtain an independent and more accurate value of Ip is underway 
[33] [47]. Once we can determine accurately the value of Ip we will be able to reject/verify 
the validity of the various cosmological models. 

Acknowledgments. We are thankful to D. Proga and P. Kernan for providing us with 
results prior to publication. 



V. FIGURE CAPTIONS 



Figure 1: Comparison between observed and calculated He I intensity ratios, on a logarithmic 
scale, for the H II region in UGC 4483 (re=16,600 K, N, ^ 100 cm'^) and the H II galaxy 
II Zw 40 (12,700 K, 225 cm~"^). The tick marks along the curves of the detailed models 
correspond to densities (from left to right) of 1, 10, 100, 10"^, 10"*, and 10^ cm~"^. The 
solid line is for the model with re=20,000 K. The plotted model of Almog & Netzer [30] 
was specially designed for comparison to Brocklehurst [44]. Note that dust extinction. Ay, 
would shift the observed line ratios towards the upper right corner and further away from 
the models. 

Figure 2: Y versus oxygen abundance (logarithmic scale relative to solar abundances) for 
the extragalactic metal-poor H II regions of Pagel et al., and the recently analyzed I Zw 18 
and UGC 4483. Such diagrams are used for extrapolation to zero metalicity (in this case, 
oxygen) to derive Ip. The recent value of Ip = 0.232 by Olive & Steigman is shown together 
with our upper bound estimated here; the same for the recent Y abundance of I Zw 18. The 
marks above the abscissa are some of the halo dwarfs from the project to obtain an accurate 
value of Ip. 
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